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A two-phase  experiment  was  conducted  on  isolated  normal  vertebral  bodies  ex- 
cised from  the  Rhesus  Macaque  at  necropsy.  The  purpose  of  this  experiment 
was  to  evaluate  the  effects  of  displacement  rate  and  vertebral  body  position 
on  mechanical  properties.  The  first  phase  of  the  experimental  tests  was 
performed  as  a factorial  experiment  in  which  12  vertebrae  from  each  of  four 
monkeys  were  loaded  and  mechanically  strained  at  each  of  three  displacement 
rates.  The  objective  of  this  phase  was  to  determine  the  effect  of  these 
factors  on  seven  mechanical  properties  and  to  determine  the  appropriate 
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regression  model  for  describing  vertebral  body  response.  In  the  second  phase 
an  additional  48  vertebrae  were  tested  at  three  different  intermediate  displace- 
ment rates  other  than  those  of  Phase  I,  again  for  the  various  vertebral  body 
positions  and  according  to  a factorial  arrangement  of  treatment  combination. 

The  results  from  Phase  II  were  compared  to  Phase  I.  The  data  from  the  two 
phases  were  combined  and  the  parameters  of  the  regression  model  were  estimated 
for  each  of  the  mechanical  properties. 


SUMMARY 


A two-phase  compression  experiment  was  conducted  on  % freshly  excised  Rhesus  monkey  vertebral  bodies. 
The  purpose  of  the  tests  was  to  evaluate  the  effects  of  displacement  rate  and  vertebral  body  position  on 
mechanical  properties.  The  objective  of  the  first  phase  was  to  determine  the  effect  of  displacement  rate  and 
vertebral  body  position  on  seven  measured  mechanical  properties.  The  appropriate  regression  model  for 
describing  the  responses  on  48  vertebrae  was  derived.  In  the  second  phase  an  additional  48  vertebrae  were 
tested  at  three  different  displacement  rates  for  vertebral  body  position  according  to  a factorial  arrangement  of 
treatment  combination.  The  results  from  Phase  II  were  compared  to  those  of  Phase  I.  The  data  from  the  two 
phases  were  combined  and  the  parameters  of  the  regression  model  were  estimated  for  each  of  the  mechanical 
properties. 
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INTRODUCTION 


The  concept  of  a dynamic  model  representing  the  axial  skeletal  system  under  the  influence  of  a short  duration, 
buttock-to-head  acceleration  is  particularly  important  when  assessing  whether  or  not  an  arbitrary  acceleration 
time  history  is  tolerable  to  man.  In  the  study  of  abrupt  accelerative  forces  and  aircrew  ejection,  attention  has 
focused  on  the  head,  spinal  column  and  pelvis.  Of  these  anatomical  components,  the  spinal  column  has 
repeatedly  been  shown  to  be  the  most  critical  to  ejection  seat  forces.  Numerous  analytical  methods  for 
describing  and  predicting  spinal  column  response  have  evolved  ranging  from  simple  lumped  parameter  single 
degree  of  freedom  to  multidegree  of  freedom,  to  continuum  models.  The  purpose  of  the  models  is  to  ascertain 
the  variation  of  risk  as  well  as  the  likelihood  of  spinal  trauma  due  to  various  acceleration  time  histories. 

Recent  advances  in  the  aerospace  sciences  have  underlined  the  fact  that  adequate  scientific  knowledge  dealing 
with  accelerative  forces,  human  spinal  response,  and  injury  potential  is  lacking.  The  shortcomings  of  current 
modeling  simulations  are  obvious.  There  exists  a need  for  more  systematic  experimental  data  collection  as  a 
basis  for  more  refined  models.  Present  day  modeling  efforts  are  far  ahead  of  the  detailed  quantitative  static  and 
dy  namic  strength  measurements  necessary  to  render  analytical  modeling  efforts  useful. 

Since  the  spinal  column  is  composed  of  discrete  components,  knowledge  of  the  biomechanical  properties  of 
individual  tissues  is  required  before  a model  can  be  represented  as  an  advanced,  discrete,  multi-parameter 
model  that  simulates  and  predicts  the  action  of  the  spine.  Data  are  required  in  the  intervertebral  disks, 
vertebral  bodies,  ligamentous  structures,  vertebral  units,  etc.  Because  of  the  relative  inaccessability  of  human 
cadaver  material  for  testing  and  the  lack  of  sufficient  quantitative  information  to  adequately  describe  the 
operational  environments,  very  limited  validation  of  human  response  and  injury  prediction  is  available  using 
current  analytical  techniques. 

Based  upon  this  information,  an  investigation  was  initiated  to  collect  static  and  dynamic  biomechanical 
strength  along  with  injury  mode  data  to  allow  the  formulation  and  experimental  corroboration  of  a subhuman 
primate  model.  Such  an  approach  would  also  authenticate  the  validation  of  current  modeling  concepts  and 
practices  and  may  identify  the  necessary  procedures  for  establishing  interspecies  scaling  relationships  that 
would  allow  extrapolation  of  kinesiologic,  kinetic,  and  injury  potential  data  between  nonhuman  primates 
and  man. 

This  investigation  deals  with  evaluating  the  effects  of  displacement  rate  and  vertebral  body  level  on  the 
mechanical  properties  of  isolated  vertebral  centra  excised  from  the  Rhesus  monkey  [ Macaca  mulatta].  The 
data  presented  herein  include  ultimate  load,  stiffness  of  deformation  to  ultimate  load  and  energy  to  ultimate 
load.  In  this  paper,  a regression  model  is  estimated  for  each  of  the  mechanical  properties. 


METHODS 

A two-phase  experiment  was  conducted  to  evaluate  the  effects  of  displacement  rate  and  position  on  the 
mechanical  properties  of  isolated  Rhesus  monkey  vertebral  bodies.  The  initial  phase  was  performed  as  a 
factorial  experiment  in  which  twelve  vertebrae  from  each  of  four  adult  Rhesus  monkeys  were  loaded  at  each  of 
three  displacement  rates.  The  objective  of  this  phase  was  to  determine  the  effect  (if  any  I of  these  factors  on 
seven  measured  mechanical  properties  and  to  determine  the  appropriate  regression  model  for  describing  and 
predicting  the  measured  responses.  In  the  second  phase,  an  additional  48  vertebrae  were  tested  at  three 
different  intermediate  displacements  other  than  those  in  Phase  I,  again  for  various  vertebral  positions 
according  to  a factorial  arrangement  of  treatment  combinations.  The  results  from  Phase  II  were  compared 
with  those  of  Phase  I.  The  data  from  the  two  phases  were  combined  and  the  parameters  of  the  regression 
model  were  estimated  for  each  of  the  mechanical  properties. 

1 

The  factorial  experiment  of  Phase  I consisted  of  compressive  tests  on  48  vertebrae  from  four  Rhesus  monkeys. 
The  twelve  vertebral  bodies  available  from  each  primate  were  grouped  into  four  sets  of  vertebral  positions  with 
three  vertebrae  in  each  set.  Position  1 (Pj  1 represents  vertebra  Tp,  T9,  and  Tip:  P2  represents T]],T]2  and 


5 


L 1 ; 1*;}  represents  L2,  L3  and  L4;  and  P4  represents  L5,  L^>  and  L7-  This  methodology  of  grouping  assumes 
that  any  differences  in  the  mechanical  properties  of  the  adjacent  vertebrae  are  negligible  in  comparison  with 
differences  between  more  distant  vertebrae.  Each  vertebral  centrum  was  subjected  to  a uniform  uniaxial 
compressive  load  to  failure  under  a constant  displacement  rate.  Within  the  respective  position  groupings,  the 
three  assigned  vertebrae  were  tested  at  one  of  three  displacement  rates:  R2  = 8.89  x 10*5  meters  per  second, 
R4  = 8.89  x 10-3  meters  pr  second,  and  R(,  = 8.89  x 10"!  meters  per  second.  IR]  and  R3  and  R5  are 
intermediate  displacement  rates  and  were  used  for  the  second  phase  of  this  experiment.  I 

The  experiment  of  Phase  II  was  conducted  similarly  to  that  of  Phase  I.  Vertebrae  from  four  different  monkeys 
were  subjected  to  compressive  loading  under  three  different  displacement  rates.  The  displacement  rates  were 
assigned  to  vertebral  position  as  In  Phase  I.  However,  the  loading  rates  in  Phase  II  were  set  at  Rj  = 8.89  x 
1 0*9  meters/ second,  R3  — 8.89  x 10-4  meters/ second  and  R5  = 8.89  x 10*2  meters/ second.  The  objective  of 
this  part  of  the  experiment  was  to  demonstrate  rates  and  to  provide  a larger  sample  size  for  estimating  the 
parameters  of  the  model. 

TEST  SPECIMEN  PREPARATION  AND  PROCEDURES 

Eight  healthy  male  Rhesus  monkeys  [ Macaca  mulatta]  were  euthanized  with  an  overdose  of  pentobarbital. 
The  subject  identification  numbers  and  body  weights  are  as  follows: 


X-96 

8.1  Kilograms 

W92 

11.4  Kilograms 

884A 

10.0  Kilograms 

34A 

8.8  Kilograms 

X90 

9.1  Kilograms 

584A 

8.2  Kilograms 

A 104 

9.3  Kilograms 

882A 

9.3  Kilograms 

Shortly  following  death,  the  vertebral  columns  were  excised  en  masse,  identified,  and  stored  in  a freezer  at 
— 30°C.  Thirty-six  hours  prior  to  testing,  the  spinal  columns  were  removed  from  the  deep  freeze  and  allowed 
to  partially  thaw.  Simultaneously,  the  individual  vertebrae  were  disarticulated  from  one  another  by  slicing 
through  the  midsection  of  the  intervertebral  disks,  the  articular  capsules  were  sectioned,  and  the  vertebral 
bodies  were  separated.  The  posterior  aspects  of  the  vertebral  bodies  were  cut  away  at  the  base  of  the  pedicles 
using  a band  saw.  (In  the  macaque,  the  transverse  processes  of  the  lumbar  and  thoracic  vertebrae  are  located 
on  the  lateral  surface  of  the  vertebral  centrum,  as  shown  in  Figure  1.)  Each  vertebral  centrum  was  cleaned  of 
all  soft  tissue  clinging  to  its  surfaces. 

All  soft  tissue  was  removed  from  the  superior  and  inferior  vertebral  body  surfaces.  The  cartilaginous  end  plate 
was  removed,  exposing  the  bony  vertebral  surface.  Both  the  superior  and  inferior  vertebral  bearing  surfaces 
were  pressed  into  an  ink  pad  and  then  onto  millimeter  paper.  Vertebral  body  bearing  area  was  determined  for 
both  inferior  and  superior  surfaces,  and  the  results  averaged.  The  length  of  the  vertebral  body  was  measured 
using  vernier  calipers. 


Figure  1.  Superior  View  of  a Vertebral  Body  from  the  Rhesus  Monkey  on  the  Left,  Human  on  the 
Right  (Note  position  of  the  transverse  process.) 
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To  promote  a uniform  load  distribution,  the  bearing  surfaces  of  each  vertebral  centrum  were  potted  in  an 
acrylic  compound.  Using  dental  acry  lic  resin  ( Acralite  88,  Keer  Manufacturing  Company,  Romulus, 

Michigan  I.  the  potting  produced  circular-shaped  pots  with  the  specimen  located  centrally.  The  diameter  of  the 
| Hit  was  subsequently  used  to  locate  the  center  of  the  specimen  coincident  with  the  loading  axis  of  the  test 
machine. 

The  vertebral  centra  were  pressed  into  the  acrylic  at  both  ends,  and  the  entire  assembly  was  placed  in  a 
V-shaped  trough  to  assure  that  both  surfaces  were  kept  parallel  and  axially  aligned  as  previously  described.**  * 
The  vertebral  centra  was  wrapped  in  a towel  soaked  • Ringer’s  solution  to  prevent  drying  while  the  acrylic 
cured.  Figure  2 shows  a vertebral  body  with  its  acrylic  bearing  surfaces  prepared  for  test. 


Figure  2.  Vertebral  Body  with  its  Acrylic  Bearing  Surfaces  Just  Prior  to  Test. 


THE  TEST  SYSTEM 

An  electrohydraulic  closed  loop  test  system  (Model  810  Material  Test  System,  MTS  Systems  Corporation, 
Minneapolis,  Minnesota  I was  used  to  strain  each  test  specimen.  The  system  is  centered  around  an  electro- 
hydraulic  closed  loop  test  machine  capable  of  being  programmed  and  controlled  in  load,  strain,  and 
displacement.  With  the  machine  in  the  displacement  control  mode,  a linear  ramp  function  was  used  to  strain 
each  test  specimen.  Testing  was  accomplished  over  five  orders  of  magnitude  of  linear  displacement  rates  18.89 
x l<r6  m/s  to  8.89  x 10"*  m/sl.  The  ultimate  ram  displacement  (specimen  deformation!  was  set  at  50%  of  the 
original  specimen  height,  up  to  a maximum  of  0.0127  meter.  The  imposed  time  dependent  displacement  and 
the  resultant  compression  loads  were  recorded.  Ram  displacements  were  measured  using  a linear  variable 
differential  transformer  ILVI^TI,  while  the  specimen  reacted  against  a four-arm  bridge  strain  gauge  load  cell. 
A multichannel  FM  magnetic  tape  recorder  and  a multichannel  transient  recorder  were  used  to  store  the  test 
results.  For  the  low-speed  tests,  load  and  displacement  data  were  recorded  directly  to  a Model  136A  Hewlett- 
Packard  recorder,  but  for  the  intermediate  and  high-speed  tests,  the  data  were  stored  in  the  digital  memory  of 
a transient  recorder  for  playback  at  reduced  speeds  into  the  X-Y  recorder.  The  FM  tape  recorder  (Electro- 
Mechanical  Research  Inc.,  Sarasota,  Florida,  Model  392AI  was  used  as  a back-up  recorder  on  the  high-speed 
tests. 


Circular  cutouts  in  the  loading  heads  of  the  test  fixture  were  utilized  to  locate  the  potted  specimens  within  the 
test  fixture.  A test  fixture  chamber  was  designed  to  contain  each  specimen  during  the  tests  while  providing  a 
steel  and  transparent  plastic  enclosure  to  be  used  for  observing  and  photographing  the  test  specimen. 

DATA  REDUCTION  METHODOLOGY 

A typical  test  curve  is  presented  in  Figure  3.  The  load  on  the  test  specimen  is  plotted  on  the  ordinate  versus  ram 
displacement  (specimen  deformation  I on  the  abscissa.  The  four  dependent  variables  were  extracted  directly 
from  the  test  curves  and  have  units  of  load  or  deformation  or  combinations  thereof.  These  variables  are 
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Figure  3.  A Typical  Load  versus  Displacement  Test  Curve. 


ultimate  load,  deformation  to  ultimate  load,  stiffness,  and  energy  to  ultimate  load.  The  first  point  on  the  load 
deformation  curve,  following  the  apparent  elastic  section,  where  the  tangent  to  the  curve  becomes  parallel  to 
the  deformation  axis  of  the  plot  (abscissa),  is  defined  as  the  ultimate  load.  The  deformation  to  this  point  is 
defined  as  the  deformation  to  idtimate  load.  The  stiffness  of  the  test  specimen  is  determined  by  fitting  the 
apparent  linear  elastic  section  of  the  load  deformation  curve  with  a least  squares  curve.  The  energy  to 
ultimate  is  defined  as  the  area  under  the  load  versus  deformation  test  curve,  from  the  point  of  zero  deformation 
up  to  the  deformation  at  ultimate  load. 

Three  additional  dependent  variables  were  determined  for  each  specimen.  These  variables  are  in  units  of 
engineering  stress  and  strain  and  are,  therefore,  dependent  upon  the  load-deformation  and  specimen  geometry 
data.  The  ultimate  engineering  stress  was  computed  from  the  ultimate  load  by  dividing  by  the  specimen 
area.  The  engineering  strain  as  ultimate  engineering  stress  was  defined  as  the  deformation  to  ultimate  load 
divided  by  the  original  length  of  the  test  specimen.  The  test  specimen  elastic  modulus  was  computed  from  the 
test  specimen  stiffness  using  the  following  relationship  with  typical  units: 

Specimen  Original  Length  (ml 

Elastic  Modulus  IPal  = Stiffness  IN /ml  9 

Specimen  Area  lm~l 


STATISTICAL  ANALYSIS 

A series  of  mechanical  compression  tests  was  conducted  on  fresh  thoracic  and  lumbar  vertebral  centra  excised 
from  the  Rhesus  monkey  [Macaca mulatto].  The  purpose  of  the  tests  was  to  identify  the  effect  of  displacement 
rate  and  vertebral  body  position  on  mechanical  strength.  In  total,  twelve  vertebral  bodies  were  tested  from 
each  vertebral  column;  of  these,  five  were  thoracic  and  seven  were  lumbar. 

The  experimental  design  was  influenced  by  several  factors.  Since  the  tests  were  destructive,  different 
displacement  rates  could  not  be  applied  to  the  same  test  specimens.  Since  the  vertebral  bodies  are  associated 
with  individual  primates  and  vertebral  body  position  is  a factor  of  interest,  the  three  displacement  rates  had  to 
be  applied  to  all  positions  in  each  primate  to  account  for  any  variation  among  primates  and  among  vertebral 
portions.  In  order  to  accommodate  the  six  displacement  rates,  it  was  assumed  that  adjacent  vertebrae  would 
exhibit  relatively  minor  variations  in  mechanical  properties  in  comparison  to  the  magnitude  of  the  effects  being 
sought.  Accordingly,  the  twelve  vertebral  centra  from  each  spinal  column  were  grouped  in  sets  of  three,  with 
each  group  representing  a position.  These  groups  were  designated  as  Pj , P2,  P3  and  P4;  where  Pi  represented 
vertebrae  Tg,  T9,  Tjo;  P2  represented  Tj  1,  T)2,  L|;  P3  represented  L2,  L3,  L4;  and  P Represented  L5, 
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Lf>  and  L7.  With  this  assumption,  each  of  the  displacement  rates  could  be  applied  to  the  subject-position 
combination.  To  the  extent  possible,  the  displacement  rates  were  balanced  across  primate  spinal  column  by 
assignments  to  particular  vertebrae  as  seen  in  the  following  matrix. 

The  spinal  columns  were  divided  into  two  groups.  Group  I consisted  of  X96  (S2I,  884A  (S3),  X90  (S4), 

A104  IS] ).  Group  II  consisted  of  W92  (Ss),  34A  (S5),  584  (S6)  and  882A  (S7).  Ri,  R2,  R3,  etc.  corresponded 
to  the  various  strain  rates  as  shown  below: 

R]  = 8.89  x 10*6  meters  per  second 
R2  = 8.89  x 10* ■>  meters  per  second 
R3  = 8.89  x 10*4  meters  per  second 
R4  = 8.89  x 10*3  meters  per  second 
R5  = 8.89  x 10*2  meters  per  second 
R(,  = 8.89  x 10*1  meters  per  second 


The  test  matrix  is  shown  in  Table  1 . For  Phase  I the  statistical  model  for  this  experiment  is  given  by: 

^ijk  — y *t*  Si  + Pj(j)  + Rfc  + SRjjj  + ejjij 

where 

Y = mechanical  properties  of  interest 
ji  = overall  mean 

Sj  = differential  effect  due  to  subject  i 
Pj(i)  = differential  effect  due  to  vertebral  position  j within  subject  i 
Rk  = differential  effect  due  to  displacement  rate  k 
SRik  = differential  joint  effect  of  subject  i and  displacement  rate  k 
eijk  — random  error  associated  with  measurement  in  subject  i,  position  j,  and  displacement  rate  k 
i = subject 
j = spinal  level 
k = strain  rate 


The  analysis  of  variance  for  the  model  is  presented  in  Table  2. 


r 


v 


TABLE  2 

GENERAL  ANALYSIS  OF  VARIANCE  TABLE  FOR  VERTEBRAL  BODY  TESTING  DATA 


Subject 

Position  (P) 

Displacement 
Rate  (R) 

RP 

Ult.  Load 

0.95 

0.95 

0.95 

— 

Deformation 

0.95 

0.95 

0.95 

— 

Stiffness 

— 

0.90 

0.95 

— 

Energy 

0.95 

0.95 

— 

— 

Ult.  Stress 

0.95 

— 

0.95 

— 

Strain  to  Ult. 

0.95 

0.95  ‘ 

0.95 

0.90 

Elastic  Mod 

0.95 

0.95 

0.95 

— 

Due  to  the  destructive  nature  of  the  tests,  replication  of  the  test  was  not  possible,  and  a pure  estimate  of  error 
variance  was  not  obtainable.  The  analysis  shown  in  the  above  table  assumed  that  the  interaction  or  joint  effect 
of  displacement  rate  at  locations  within  subject  was  negligible,  thus  permitting  this  test  to  be  used  as  the 
estimate  of  error  variance.  (The  plausibility  of  the  assumption  was  verified  by  the  data  of  this  experiment.  I To 
calculate  the  mean  squares  of  Table  3,  the  test  specimens  were  considered  as  being  random  samples  from  a 
population  where  position  and  displacement  rate  were  considered  as  fixed  effects. 

To  test  for  a significant  difference  between  differential  effects  of  the  independent  variables  in  an  analysis  of 
variance,  the  test  statistically  defined  by  the  appropriate  F ratio  of  observed  mean  squares  was  calculated. 
Under  the  null  hypothesis  of  no  difference,  the  variance  term  (a^)in  the  expected  mean  squares  of  that  effect 
is  zero,  and  the  expected  F ratio  is  unity.  Large  F ratios  calculated  from  observed  data  would  indicate  that  the 
effect  is  present  where  the  largeness  can  be  determined  for  a given  level  of  confidence  from  a table  of  critical  F 
values.  The  critical  F value  was  a function  of  the  degrees  of  freedom  in  the  numerator  and  the  denominator  of 
the  F ratio.  Note  that  for  this  experiment,  the  effects  of  subjects,  positions  and  the  subject-by-displacement- 
rate  interaction  were  determined  by  making  a ratio  of  their  observed  mean  square  with  that  of  the  error  term. 
The  displacement  rate  effect,  however,  was  determined  by  making  a ratio  of  its  observed  mean  square  with 
that  of  the  subject  by  displacement  rate  interaction.  Once  degrees  of  freedom  are  analogous  to  sample  size,  it 
can  be  seen  that  a smaller  sample  size  is  available  for  testing  for  displacement  rate  effect  than  for  testing  the 
effect  of  the  other  sources  of  variation.  Increasing  the  number  of  subjects  would  increase  the  sample  size. 

The  experiment  of  Phase  II  was  conducted  similarly  to  that  of  Phase  I.  Vertebral  centra  from  four  different 
monkeys  were  subjected  to  compressive  loading  under  three  different  displacement  rates.  The  loading  rates  in 
this  Phase  were  set  at  Rl  = 8.9  x 10'^  m/s,  R3  = 8.9  x 10"^  m/s,  and  R5  = 8.9  x 10"^  m/s.  The  objective  of 
this  experiment  was  to  demonstrate  that  the  form  of  the  model  derived  in  Phase  I would  apply  to  other 
displacement  rates  and  to  provide  a larger  sample  size  for  estimating  the  parameters  of  the  model. 

The  results  from  one  of  the  test  subjects  W92  (Ss)  was  significantly  different  from  the  oth-  r seven  subjects 
tested  and  as  a result  were  not  included  in  this  analysis. 

The  regression  analysis  of  the  Phase  I data  indicated  that  an  appropriate  model  for  predicting  any  of  the 
mechanical  responses  in  Phase  II  has  the  form: 

Yij  = ap;  + bx  + e;j 

where 
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Yjj  = jth  response  in  position  grouping  i 
Bp,  = intercept  for  position  grouping  i 
x = log  ip  of  displacement  rate 
b = response  slope 
ejj  = random  error 

Note  that  the  qualitative  variable  of  vertebral  position  is  modeled  by  generating  equations  for  different 
positions.  The  Phase  I statistical  analysis  indicated  that  the  effect  due  to  position  was  not  significant; 
therefore,  only  a single  equation  was  derived  for  positions  within  the  vertebral  column.  Similarly,  when 
subgroupings  of  position  displayed  significant  differences  between  subgroups  but  the  difference  within  a 
subgrouping  was  not  significant,  a different  constant  was  calculated  for  each  subgroup. 

TABLE  3 

PHASE  I - ULTIMATE  LOAD  (N) 

AND  ANALYSIS  OF  VARIANCE  TABLE 


P1 

P2 

P3 

P4 

H2 

S1 

2447 

4092 

4849 

5516 

S2 

2113 

2936 

4048 

4448 

S3 

2157 

3292 

4715 

5204 

S4 

2068 

2580 

2936 

3203 

R4 

S1 

3336 

3870 

5783 

6494 

S2 

2891 

* 

4893 

4671 

5694 

S3 

3190 

4444 

5738 

6583 

S4 

2091 

2313 

3205 

4003 

R6 

S1 

4448 

5961 

7206 

8096 

S2 

3692 

6032 

5516 

8452 

S3 

2936 

4003 

* 

6672 

8185 

S. 

4 

2802 

3718 

4359 

5738 

‘Estimated  value  of  missing  data  point. 


Analysis  of  Variance  Table 


Source  of 
Variation 

Degrees  of 
Freedom 

Mean 

Square 

F 

Ratio 

Level  of 
Confidence 

3 

8.478  x 10* 

23.24 

3 

15.438  x 10® 

42.26 

BmI 

2 

21.298  x 10® 

58.31 

0.95 

6 

0.612  x 10® 

1.67 



31 

0.365  x 10 
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The  ultimate  loads  for  Phase  II  are  shown  in  Table  4.  The  effects  due  to  both  displacement  rate  and  vertebral 
body  position  were  also  found  to  be  significant.  The  ultimate  load  for  R5  was  greater  than  R3,  which  was  in 
turn  greater  than  R 1 . 


TABLE  4 

PHASE  II  RESULTS  - ULTIMATE  LOADS  (N) 


P2 

p3 

P4 

R1 

S5 

2157 

2713 

3781 

3781 

s6 

2269 

2758 

3136 

3247 

S7 

2891 

3514 

3781 

5471 

S8 

867 

1090 

1601 

1668 

R3 

S5 

2891 

3737 

4537 

4671 

s6 

2313 

2713 

3781 

4181 

s7 

4300 

3914 

5605 

6939 

S8 

1223 

1334 

2535 

1824 

R5 

S5 

1935 

3648 

5071 

mm 

S6 

2580 

3959 

5115 

ESI 

s7 

4212 

5783 

6850 

8007 

S8 

1646 

1913 

2269 

2402 
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RESULTS  AND  DISCUSSIONS 


The  values  of  the  mechanical  properties  that  were  reduced  Irom  the  test  data  are  given  in  Appendix  I and 
Appendix  II. 

ULTIMATE  LOAD 

The  ultimate  loads  along  with  the  calculated  analysis  of  variance  are  shown  in  Table  3.  The  effects  due  to  both 
displacement  rate  and  vertebral  position  were  found  to  be  significant.  In  particular,  the  average  ultimate  load 
for  R(>  was  significantly  greater  than  that  for  R4,  which  in  turn  was  significantly  greater  than  the  average 
ultimate  load  for  R2-  Similarly,  the  average  ultimate  loads  for  the  four  positions  were  all  significantly  different 
with  P4>P3>P2 >Pl . Figure  4 identifies  average  load  as  a function  of  displacement  rate  for  each  position.  The 
results  indicate  that  ultimate  load  is  a linear  function  of  log  displacement  rates,  with  different  constants  being  a 
function  of  vertebral  level  or  position. 


P.  P*  P*  P4 

POSITION 


Figure  4.  Average  Ultimate  Load  for  Position  by  Displacement  Rate  Combinations. 
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DEFORMATION  TO  ULTIMATE  LOAD 

The  values  of  deformation  to  ultimate  load  observed  during  Phase  I and  their  analysis  of  variance  summary 
are  presented  in  Table  5.  The  effects  due  to  displacement  rate  and  position  were  significant.  The  average 
deformation  for  Rft  was  significantly  greater  than  that  of  R4,  which  in  turn  was  significantly  greater  than  that 
of  R2-  The  average  deformation  for  Pj  was  significantly  less  than  the  other  three,  and  that  of  P4  was 
significantly  greater  than  the  others.  No  significant  difference  was  noted  between  P2  and  P3.  The  results 
shown  in  Figure  5 represent  the  average  deformation  for  combination  of  position  and  displacement  rate.  The 
linear  change  of  deformation  with  log  displacement  rate  is  not  as  clearly  defined  as  in  the  ultimate  load  values. 
The  data  show  considerably  more  scatter. 


P, 


P* 

POSITION 


PS 


P4 


Figure  5.  Average  Deformation  to  Ultimate  Load  for  Position  by  Displacement  Rate  Combinations. 
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TABLE  5 


PHASE  I RESULTS  - DEFORMATION  TO  ULTIMATE  LOAD  (m  x 10'3) 


Pi 

p2 

p3 

P4 

*2 

Si 

4 . 890 

4.648 

5.080 

6.096 

s2 

2.699 

2.540 

3.937 

4.318 

s3 

2.540 

2.350 

2.667 

3.175 

s4 

4.445 

2.883 

2.667 

4.763 

R4 

S1 

4.763 

4.604 

4.445 

5.080 

s2 

3.016 

3.  366 

2.223 

2.667 

s3 

1.868* 

2.223 

3.493 

s4 

1.270 

3.810 

3.810 

r6 

S1 

1.778 

4.128 

4.445 

5.715 

s2 

2.477 

3.175 

2.527 

5.398 

S3 

2.032 

1.937 

2.527 

2.413 

s4 

2.381 

2.767* 

2.032 

3.810 

‘Estimated  value  of  missing  data  point. 
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Analysis  of  Variance  Table 


Source  of 
Variation 

Degrees  of 
Freedom 

HESTE 

m 

S 

3 

9.690 

18.94 

0.95 

R 

2 

1.693 

3.31 

0.95 

P 

3 

4.253 

8.31 

0.95 

RP 

6 

0.671 

1.31 

— 

Error 

31 

0.512 
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The  Phase  II  results  with  respect  to  deformation  to  ultimate  load  are  similar  to  that  of  the  above;  the  average 
deformation  for  R",  was  significantly  greater  than  that  of  R3,  which  in  turn  was  conclusively  above  that  of  R ] . 
These  results  are  shown  in  Table  6. 


TABLE  6 

PHASE  II  RESULTS  - DEFORMATION  TO  ULTIMATE  LOAD  (m  x HP3) 


P1 

p2 

P3 

, 

P4 

m 

s5 

5.080 

4.763 

6.350 

6.223 

■ 

s6 

1.778 

2.159 

3.810 

8.128 

s7 

4.445 

5.080 

3.366 

6.985 

■ 

S8 

2.159 

2.540 

4.572 

1.778 

R3 

S5 

5.398 

4.76  3 

8.255 

6.668 

S6 

2.096 

2.858 

3.493 

4.445 

s7 

— 

3.175 

5.715 

5.906 

S8 

2.191 

3.239 

4.445 

1.905 

R5 

s5 

1.429 

4.604 

5.715 

4.572 

S6 

1.588 

2.699 

3.493 

5.080 

S7 

— 

5.080 

3.937 

5.080 

S8 

3.651 

2.858 

1.969 

1.588 

STIFFNESS 

The  stiffness  values  and  their  analysis  of  variance  are  shown  in  Table  7.  Displacement  rate  was  significant  with 
higher  rates  of  loading  resulting  in  larger  stiffness  values.  The  position  effect  was  significant  at  the  90  percent 
confidence  level,  however  not  at  the  95  percent  level  of  confidence.  Since  the  means  for  P]  and  P2  (2.62  and 
2.81  N/m'x  l(/\  respectively!  were  approximately  equal,  separate  regression  models  were  derived  for  each  of 
the  two  position  groupings.  The  results  are  presented  in  Figure  6,  which  provides  average  stiffness  for 
combinations  of  displacement  rates  and  position.  The  response  of  stiffness  due  to  the  controlled  i actors  is  not 
always  consistent  for  the  individual  position  grouping  or  displacemeni  • ates,  but  the  significant  trends  can  be 
discerned.  A similar  trend  was  noted  to  be  present  in  the  Phase  II  stiffness  data  shown  in  Table  8. 


17 


TABLE  7 


PHASE  I RESULTS  - STIFFNESS  (N/mxlO6) 


Pi 

P2 

P3 

P4 

*2 

si 

1.067 

1.879 

1.728 

2.179 

S2 

0.885 

1.492 

1.788 

1.230 

S3 

1.023 

1.804 

2.257 

2.058 

S4 

6.273 

1.057 

1.696 

1.584 

m 

S1 

1.268 

1.463 

3.395 

3 . 1 5 i 

s2 

1.226 

1.783 

4.041 

4.028 

1 

s3 

1.809* 

1.611 

3.573 

3.065 

■ 

s4 

2.038 

0.866 

1.051 

2.335 

Si 

2.709 

3.379 

3.940 

s2 

1.121 

2.001 

5.004 

S 3 

1.821 

4.641 

3.244 

S4 

2.371 

3.337* 



3.503 

*Estimated  value  of  missing  data  point. 


Analysis  of  Variance  Table 


Source  of 

Variation 

Degrees  of 

Freedom 

Mean 

Square 

F 

Ratio 

Level  of 
Confidence 

S 

3 

0.448 

0.34 



R 

3 

8.115 

6.15 

0.95 

P 

2 

3.126 

2.37 

0.90 

RP 

6 

1.831 

1.39 

— 

Error 

31 

1.  319 
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STIFFNESS  (N/m  x 10 


TABLE  8 


PHASE  II  RESULTS  - STIFFNESS  (N/mx  106) 


P1 

P2 

— ; 

p3 

P4 

R1 

S5 

0.711 

0.861 

1.868 

1.373 

S6 

1.576 

1.560 

1.101 

0.881 

S7 

1.061 

0.799 

1.653 

1.217 

s8 

0.458 

0.494 

0.457 

1.066 

R3 

s5 

0.943 

1.907 

1.284 

1.541 

S6 

1.261 

1.051 

1.423 

1.479 

S7 

— 

1.910 

1.615 

2.215 

s8 

0.631 

0.480 

0.978 

1.194 

r5 

S5 

1.541 

1.611 

1.500 

S6 

1.821 

1.681 

"1 

1.691 

s7 

— 

1.416 

3.121 

2.802 

S8 



0.490 

1.090 

1.168 

1.518 
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ENERGY  TO  ULTIMATE  LOAD 

The  values  of  energy  to  ultimate  load  obtained  during  the  Phase  I tests  are  shown  in  Table  9 along  with  the 
analysis  of  variance  summary.  The  effect  due  to  displacement  rate  was  not  significant,  but  the  effect  due  to 
vertebral  centrum  position  within  the  spinal  column  was  significant.  These  results  are  discernible  from 
Figure  7, which  presents  average  energy  for  combination  of  displacement  rate  at  position.  Like  results  are  noted 
in  the  Phase  II  energy  to  ultimate  load  data  shown  in  Table  10. 


2 


J 1 L_ 

Rj 

DISPLACEMENT  RATE 


P* 


POSITION 


P4 


Figure  7.  Average  Energy  to  Ultimate  Load  for  Position  by  Displacement  Rate  Combinations. 
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TABLE  9 


PHASE  I RESULTS  - ENERGY  TO  ULTIMATE  LOAD  (J) 


pl 

P2 

p3 

P4 

r2 

Si 

7.987 

24.049 

s2 

2.881 

3.994 

3.717 

9.869 

S 3 

2.870 

4.429 

6.463 

8.920 

S4 

5.357 

3.850 

10.739 

r4 

S1 

9.736 

10.996 

21.942 

s2 

4.316 

9.276 

6.999 

10.344 

S3 

1.629* 

6.022 

7.338 

14.349 

S4 

1.370 

3.265 

2.146 

10.434 

R6 

S1 

3.305 

16.597 

22.134 

27.580 

s2 

3.839 

7.355 

3.921 

30.834 

S3 

3. 101 

5.268 

7.380 

12.248 

— 

s4 

4.384 

— 

5.308* 

5.259 

15.389 

‘Estimated  value  of  missing  data  point. 


Analysis  of  Variance  Table 


WnSBSBSM 

mSSSSSm 

mmssM 

F 

Ratio 

Level  of 
Confidence 

s 

3 

386.67 

29.15 

0.95 

R 

2 

29.26 

2.2] 

— 

P 

3 

316.18 

23.84 

0.95 

RP 

6 

18.92 

1.43 

— 

Error 

31 

13.26 
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TABLE  10 


PHASE  II  RESULTS  - ENERGY  TO  ULTIMATE  LOAD  (J) 


P1 

P2 

P3 

P4 

-x  1 

S5 

6.799 

7.819 

17.451 

16.468 

S6 

2.642 

2.949 

6.672 

10.214 

s7 

7.453 

7.666 

7.700 

22.563 

s8 

0.859 

1.408 

4.000 

1.347 

R3 

S5 

9.827 

12.835 

25.970 

21.173 

S6 

2.583 

3.632 

6.728 

11.988 

S? 

— 

7.604 

19.315 

26.269 

s8 

1.637 

2.110 

7.355 

1.638 

R5 

S5 

1.212 

11.462 

18.959 

12.530 

s6 

1.840 

4.638 

9.756 

16.965 

S7 

— 

17.982 

16.383 

25.681 

S8 

2.781 

3.197 

1.667 

1.678 
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ULTIMATE  ENGINEERING  STRESS 

The  values  of  ultimate  engineering  stress  and  their  analysis  of  variance  summary  are  presented  in  Table  1 1 . 
The  effect  of  displacement  rate  is  significant  with  higher  values  of  ultimate  stress  resulting  from  the  faster 
displacement  rates.  The  position  effect,  however,  was  not  significant.  These  results  can  be  seen  in  Figure  8, 
which  presents  average  ultimate  stress  for  combinations  of  displacement  rate  and  position.  Phase  II  results  are 
presented  in  Table  12.  The  results  are  similar  to  the  above  findings. 


TABLE  11 

PHASE  I RESULTS  - ULTIMATE  ENGINEERING  STRESS  (Pax  10*) 


P1 

P2 

p3 

p4 

r2 

19.056 

19.049 

19.987 

21.644 

14.621 

14.220 

15.454 

15.494 

14.995 

18.621 

18.456 

17.966 

11.963 

9.874 

9.230 

10.257 

r4 

S1 

21.457 

21.577 

25.177 

23.410 

s2 

17.998 

19.497 

17.237 

18.048 

S3 

22.300* 

21.465 

22.517 

21.437 

s4 

15.579 

10.545 

11.750 

12.215 

21.669 

29.707 

29.163 

30.756 

20.810 

22.638 

23.686 

25.486 

26.768 

26.294 

29.133 

29.850 

18.405 

17.876* 

14.531 

16.941 

‘Estimated  value  of  missing  data  point. 


Analysis  of  Variance  Table 


■iSBM 

wmnttrjwm 

M 

Level  of 
Confidence 

S 

3 

259.07 

75.24 

0.95 

R 

2 

280.41 

81.44 

0.95 

mm 

3 

5.67 

1.65 

— 

6 

3.56 

1.03 

— 

31 

3.44 
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Figure  8.  Average  Ultimate  Engineering  Stress  for  Position  by  Displacement  Rate  Combinations. 


TABLE  12 


PHASE  II  RESULTS  - ULTIMATE  ENGINEERING  STRESS  (Pa  x 10®) 


pl 

P2 

P3 

P4 

R1 

s5 

15.699 

13.835 

15.670 

13.442 

S6 

15.288 

12.573 

12.368 

12.939 

S7 

15.454 

14.801 

13.790 

16.793 

S8 

7.307 

5.312 

6.159 

5.409 

R3 

S5 

16.912 

17.288 

15.663 

16.681 

S6 

16.446 

15.021 

15.066 

15.693 

s7 

— 

20.292 

18.061 

20.332 

s8 

7.645 

6.249 

9.380 

6.132 

R5 

S5 

15.703 

' 

20.559 

21.358 

17.966 

S6 

18.774 

20.731 

19.922 

20.520 

S7 

— 

22.692 

22.591 

23.068 

s8 

9.277 

7.741 

7.745 

7.708 

ENGINEERING  STRAIN  TO  ULTIMATE  STRESS 

The  strain  values  obtained  during  the  Phase  I tests  are  presented  in  Table  13  along  with  their  analysis  of 
variance.  The  displacement  rate  and  position  effects  are  significant  at  the  95%  level  of  confidence,  while  the 
joint  effect  of  displacement  rate  at  position  is  significant  at  the  90%  level  of  confidence.  Since  this  is  the  only 
response  for  which  the  interaction  is  significant  and  since  it  is  not  significant  at  a high  level  of  confidence,  the 
interaction  term  was  removed  in  modeling  the  strain  response  as  a function  of  displacement  rate  at  position. 
The  average  strain  measurements  did  not  display  a consistent  pattern;  this  is  shown  in  Figure  9.  There  is  a 
decreasing  average  strain  for  increasing  displacement  rate  and  a trend  toward  lower  strains  for  P3  and  P4  than 
for  Pi  and  P2.  As  a result,  strain  was  modeled  in  terms  of  log  displacement  rate  with  Pi  and  P2  in  one 
subgrouping  and  P3  and  P4  in  a record  subgrouping.  Phase  II  results  as  shown  in  Table  14  reflected  similar 
findings. 
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TABLE  13 


PHASE  I RESULTS  - ENGINEERING  STRAIN  TO  ULTIMATE  STRESS  (m/ml 


Pi 

P2 

— 

P3 

P4 

m 

S1 

0.507 

0.306 

H 

s2 

0.280 

0.182 

s3 

0.265 

0.182 

■ 

s4 

0.378 

0.185 

m 

0.469 

0.376 

0.266 

0.300 

ic 

0.282 

0.228 

0.110 

0.124 

0.229* 

0.203 

0.111 

0.168 

■ 

0.130 

0.238 

0.206 

0.179 

R6 

S1 

0.161 

0.308 

0.232 

0.293 

s2 

0.191 

0.191 

0.122 

0.269 

s3 

0.218 

0.159 

0.153 

0.112 

s4 

0.233 

0.188* 

0.104 

0.176 

‘Estimated  value  of  missing  data  point. 


Analysis  of  Variance  Table 


wnT.  "232 

wsm 

WBSEEEM 

S 

3 

0.05007 

17.48 

0.95 

R 

2 

0.01131 

3.95 

0.95 

P 

3 

0.02294 

8.01 

0.95 

RP 

6 

0.00638 

2.23 

0.90 

Error 

31 

0.00286 
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R2  r4  r6  pl  p2  p3  p4 

DISPLACEMENT  RATE  POSITION 


Figure  9.  Average  Engineering  Strain  to  Ultimate  for  Position  by  Displacement  Rate  Combinations. 
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TABLE  14 


PHASE  II  RESULTS  - ENGINEERING  STRAIN  TO  ULTIMATE  STRESS  (m/m) 


ELASTIC  MODULUS  IPo  x 10 


ELASTIC  MODULUS 

The  Phase  I values  of  elastic  modulus  and  their  analysis  of  variance  are  presented  in  Table  15.  Both  the 
displacement  rates  and  position  effects  are  significant.  The  elastic  modulus  increases  with  increasing 
displacement  rate  as  is  apparent  from  Figure  10.  The  increase  is  approximately  linear  with  log  displacement 
rate.  The  average  elastic  modulus  for  Pi  is  significantly  less  than  that  of  P2,  which  in  turn  is  significantly  less 
than  that  of  P3  and  P4.  The  differences  in  average  elastic  modulus  between  P3  and  P4  is  not  significant. 
These  trends  are  readily  observable  in  Figure  10  and  are  also  similar  as  noted  in  Table  16,  which  presents 
Phase  II  results. 


DISPLACEMENT  RATE  POSITION 


Figure  10.  Average  Elastic  Modulus  for  Position  by  Displacement  Rate  Combinations. 
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TABLE  IS 


PHASE  I RESULTS  - ELASTIC  MODULUS  (Pa  x 10®) 


R1 

r2 

R3 

P4 

r2 

S1 

80.2 

133.1 

123.2 

172.7 

s2 

58.9 

100.6 

127.1 

93.1 

S3 

68.3 

132.0 

162.9 

156.9 

S4 

42.7 

63.2 

107.8 

99.3 

R4 

S1 

82.8 

99.9 

247.5 

192.4 

s2 

81.6 

104.8 

300 . 3 

274.8 

s3 

157.1* 

115.7 

279.9 

207.1 

S4 

148.1 

50.0 

75.5 

R6 

Si 

188.9 

225.8 

305.7 

s2 

82.0 

125.0 

377.3 

s3 

154.4 

370.4 

233. 8 

■ 

S4 

138.1 

188.3* 

228.3 

‘Estimated  value  of  missing  data  point. 
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Analysis  of  Variance  Table 


Degrees  of 
Freedom 

Mean 

Square 

F 

Ratio 

Level  of 
Confidence 

■Hi 

3 

12498 

3.47 

0.95 

2 

59861 

16.63 

0.95 

HLH 

3 

30158 

8.38 

0.95 

6 

3717 

1.03 

— 

31 

3599 

31 


TABLE  16 


PHASE  II  RESULTS  - ELASTIC  MODULUS  (PaxlO6) 


P1 

P2 

p3 

P4 

r 

Rl 

S5 

56.7 

72.0 

162.4 

111.0 

s6 

129.0 

117.6 

89.1 

62.9 

S7 

69.9 

52.3 

114.8 

85.8 

s8 

36.4 

29.2 

31.1 

70.4 

r3 

s5 

70.3 

160.3 

108.8 

90.5 

S6 

88.8 

78.2 

102.9 

121.4 

s7 

— 

135.8 

210.2 

147.5 

S8 

38.9 

32.0 

67.1 

83.6 

R5 

s5 

128.0 

130.6 

133.5 

128.4 

SS 

140.4 

130 . 8 

137.9 

128.8 

S7 

— 

97.3 

110.4 

156.7 

S8 

28.8 

68.3 

89.9 

94.3 

The  mechanical  property  responses  of  interest  just  described  were:  ultimate  load,  deformation  to  ultimate 
load,  stiffness,  energy  to  ultimate  load,  ultimate  engineering  stress,  engineering  strain  to  ultimate  stress,  and 
elastic  modulus.  For  each  of  these  dependent  variables,  an  analysis  of  variance  was  performed  on  the  resulting 
measurements.  Table  1 7 summarizes  the  results  of  these  analyses.  Since  variation  between  monkeys  was 
expected,  this  source  of  variation  was  isolated  in  the  analyses  to  provide  better  discrimination  among  the  levels 
of  the  other  factors.  A significant  effect  for  subject  (monkeys)  indicates  that  differences  in  response  were 
observed  for  the  monkeys,  but  this  observation  has  no  other  practical  significance. 
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TABLE  17 


SUMMARY  OF  REGRESSION  RESULTS  FOR  RHESUS  MONKEYS 


Physical 

Property 

R 

S(e) 

Slope 

(b) 

Position 

Grouping 

3Pi 

Ultimate 

Load  (N) 

0.83 

936 

448 

P1 

3 

3014 

4230 

5132 

5965 

Deformation 
(m  x 10~3) 

0.55 

1.317 

-0.271 

pi 

p2-p3 

p4 

2.796 

3.533 

4.754 

Stiffness 
(N/m  x 106) 

0.56 

1.953 

0.320 

P1  'P2 
p3'p4 

1.978 

2.613 

Energy (J) 

0.65 

5.64 

0.399 

pl 

p2 

ij 

4.836 

8.424 

11.667 

17.449 

Ultimate 

Eng.  Stress 
(Pa  x 106) 

0.68 

3.657 

1.911 

pl'p2'p3' 

p4 

19.911 

Strain 

(m/m) 

0.41 

0.086 

-0.0159 

pl'p2 

P3'P4 

0.247 

0.212 

Elastic 

Modulus 

0.68 

58.5 

24.7 

pl 

p2 

P3'P4 

119.6 

145.1 

188.2 

The  entries  in  Table  1 7 represent  the  level  of  confidence  of  rejecting  the  hypothesis  of  no  effect  due  to  the 
factor . Levels  of  confidence  greater  than  0.95  were  not  considered.  Levels  of  confidence  as  low  as  0.90  were 
reported  since  consistent  trends  were  noted  in  the  data  and  the  large  degree  of  scatter  could  prevent  the  trends 
from  being  statistically  significant  at  a higher  level  of  confidence.  The  absence  of  an  entry  indicates  the  effect 
is  not  significant  at  the  90%  level  of  confidence. 

THE  REGRESSION  MODEL 

The  qualitative  variable  of  vertebral  position  is  modeled  by  generating  equations  for  different  positions. 
However,  when  the  Phase  I analysis  indicated  that  the  effect  due  to  position  was  not  significant,  only  one 
equation  was  derived  for  all  positions.  Similarly,  when  subgroupings  of  position  displayed  significant 
differences  between  subgroups  but  the  differences  within  a subgroup  were  not  significant,  a different  constant 
was  calculated  for  each  subgroup. 
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The  above  model  was  fit  to  the  combined  data  sets  of  Phase  I and  II,  and  the  results  are  summarized  in  Table 
17.  For  each  of  the  physical  properties,  the  table  contains  a single  value  for  the  correlation  coefficient,  the 
standard  error,  and  the  slope  (coefficient  of  log  displacement  rate).  The  groupings  of  vertebral  position  are  also 
indicated  along  with  the  intercept  for  each  group.  The  percentage  of  total  variation  of  the  response  mechanical 
properties  explained  by  regression  (r^)  varied  from  a high  of  69%  for  ultimate  load  down  to  a low  of  17%  for 
engineering  strain  to  ultimate  stress.  Thus,  although  the  regression  equations  do  provide  a significant 
correlation  with  the  measured  mechanical  properties,  considerable  unexplained  variability  due  to  uncontrolled 
and/or  unmeasured  factors  still  exists. 

Figures  1 1 through  17  display  the  regression  results  for  the  first  listed  position  grouping  of  Table  17.  Load, 
deformation,  energy,  and  elastic  modulus  are  presented  for  Pj  (vertebral  bodies  Tg,  T9  and  T]()|.  Stiffness 
and  strain  are  presented  for  Pi  and  P2  (vertebral  bodies  Tg,  T9,  T10,  Tn,  T12  and  Ll  I since  no  significant 
difference  in  these  responses  was  detected  for  these  position  groupings.  The  figure  for  stress  (Figure  1 1 1 is 
applicable  to  all  positions  since  this  factor  did  not  significantly  influence  the  ultimate  stress  determinations. 
Since  the  only  difference  between  the  regression  equations  for  different  positions  is  an  additive  constant,  the 
presented  curves  are  representative  of  all  positions. 

The  five  curves  in  each  of  these  figures  are  estimates  of  percentiles  of  the  individual  responses  at  a particular 
displacement.  That  is,  25%  of  the  response  values  from  individual  vertebrae  would  be  expected  to  fall  below  the 
25th  percentile  line.  Further,  90%  of  individual  responses  would  fall  between  the  5th  and  95th  percentile  lines. 
These  percentile  lines  (or  others  if  desired)  are  based  on  the  assumption  of  the  random  errors  (ejj)  being 
normally  distributed  with  equal  variance  for  all  displacement  rates  and  can  be  generated  from  the  data  of 
Table  17. 

Figures  1 1 through  17  also  present  the  individual  observed  test  values  after  correcting  for  applicable  position 
effects.  The  correction  is  an  additive  constant  equal  to  the  difference  in  ap;  values  from  Table  17.  Two 
conclusions  can  be  drawn  from  an  examination  of  the  plots  of  the  data.  First,  the  two  phases  of  the  experiment 
resulted  in  a consistent  relationship  between  the  response  properties  over  the  range  of  displacement  rates.  The 
results  of  the  two  phases  were  in  agreement.  Second,  the  statistical  assumptions  required  for  the  regression 
analysis  are  not  grossly  violated,  and  the  percentile  lines  contain  approximately  the  correct  percentages  of 
data  points. 
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LOG  R ( R in  m/  min  ) 


Figure  11.  Regression  Results  for  Stress  vs.  Displacement  Rate  — Pi. 


LOG  R (R  in  m/min) 


Figure  12.  Regression  Results  for  Deformation  vs.  Displacement  Rate  — Pj. 


LOG  R (R  in  m/min) 


Figure  13.  Regression  Results  for  Stiffness  vs.  Displacement  Rate  — (Pj  and  P2I 
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Figure  15.  Regression  Results  for  Stress  vs.  Displacement  Rate  — All  positions. 
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LOG  R (R  in  m/min) 


Figure  16.  Regression  Results  for  Strain  vs.  Displacement  Rate  — Pj  and  P2- 
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T — V 


Figure  17.  Regression  Results  for  Elastic  Modulus  vs.  Displacement  Rate  — P|. 
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APPENDIX  I 


ENGINEERING  DATA  OBTAINED  FROM  RHESUS 
MONKEY  SINGLE  VERTEBRAL  BODIES 


Speclman 

No. 

VB 

No. 

Displace. 

Rate 

(m/min) 

Ult. 

Load 

(N) 

Def . to 

Ult .Load 
(MxlO-3) 

Stiffness 

(N/mxl06) 

Energy  to 
Ult.  Load 
(J) 

884A 

T8 

184 

53.34 

2936 

2.032 

1.821 

3.101 

884a 

T9 

185 

5.334xl0-3 

2157 

2157 

2.540 

2.870 

884A 

T10 

186 

0.5334 

— 

— 

— 

— 

884A 

Til 

187 

53.34 

4003 

1.937 

4.641 

5.268 

884A 

T12 

188 

5.334x10-3 

3292 

2.350 

1.804 

4.429 

884A 

LI 

189 

0.5334 

4444 

2.985 

1.611 

6.022 

884A 

L2 

190 

53.34 

6672 

2.527 

3.244 

7.830 

884A 

L3 

191 

5.334x10-3 

4715 

2.667 

2.257 

6.463 

884A 

L4 

192 

0.5334 

5738 

2.223 

3.573 

7.338 

884A 

L5 

193 

53.34 

8185 

2.413 

5.511 

12.248 

884A 

L6 

194 

5.334xl0-3 

5204 

3.175 

2.058 

8.920 

884A 

L7 

195 

0.5334 

6583 

3.493 

3.065 

14.349 

A104 

T8 

196 

5.334x10-3 

2447 

4.890 

1.067 

7.987 

A104 

T9 

197 

0.5334 

3336 

4.763 

1.268 

9.736 

A104 

T10 

198 

53.34 

4448 

1.778 

2.709 

3.305 

A104 

Til 

199 

0.5334 

3870 

4.604 

1.463 

10.996 

A104 

T12 

200 

53.34 

5961 

4.128 

3.379 

16.597 

A104 

LI 

201 

5.334x10-3 

4092 

4.648 

1.879 

14.259 

A104 

L2 

202 

0.5334 

5783 

4.445 

3.395 

18.241 

A104 

L3 

203 

5.334x10-3 

4849 

5.080 

1.728 

16.456 

A104 

L4 

204 

53.34 

7206 

4.445 

3.940 

22.134 

A104 

L5 

205 

5. 334xl0-3 

5516 

6.096 

2.179 

24.049 

A104 

L6 

206 

53.34 

8096 

5.715 

2.335 

27.580 

A104 

L7 

207 

0.5334 

6494 

5.080 

3.152 

21.942 

X90 

T8 

208 

0.5334 

2091 

1.270 

2.038 

1.370 

X90 

T9 

209 

53.34 

2802 

2.381 

2.371 

4.384 

X90 

T10 

210 

5. 334x10" i 

2068 

4.445 

6.273 

5.357 

X90 

Til 

211 

0.5334 

2313 

3.016 

0.866 

3.265 

X90 

T12 

212 

53.34 

— 

— 

— 

— 

X90 

LI 

213 

5.334x10-3 

2580 

2.883 

1.057 

3.850 

X90 

L2 

214 

0.5334 

3025 

3.810 

1.051 

2.146 

X90 

L3 

215 

53.34 

4359 

2.032 

3.503 

5.259 

X90 

L4 

216 

5.334xl0-3 

2936 

2.667 

1.696 

4.678 

X90 

L5 

217 

0.5334 

4003 

3.810 

2.335 

10.434 

X90 

L6 

218 

53.34 

5738 

3.810 

3.585 

15.389 

X90 

L7 

219 

5.334x10-3 

3203 

4.763 

1.584 

10.739 

X96 

T8 

220 

5.334x10-3 

2113 

2.699 

0.885 

2.881 

X96 

T9 

221 

0.5334 

2891 

3.016 

1.226 

4.316 

X96 

T10 

222 

53.34 

3692 

2.477 

1.121 

3.839 

X96 

Til 

223 

5.334x10-3 

2936 

2.540 

1.492 

3.994 

X96 

T12 

224 

0.5334 

4893 

3.366 

1.783 

9.276 

X96 

LI 

225 

53.34 

6032 

3.175 

2.001 

7.355 

X96 

L2 

226 

5.334x10-3 

4048 

3.937 

1.788 

8.717 

APPENDIX  I (cont) 


ENGINEERING  DATA  OBTAINED  FROM  RHESUS 
MONKEY  SINGLE  VERTEBRAL  BODIES 


Specimen 

No. 

VB 

No. 

Displace. 

Rate 

(m/min) 

Ult. 

Load 

(N) 

Def . to 

Ult .Load 
(mxl0-3) 

Stiffness 

(N/mxl0^) 

Energy  to 
Ult.  Load 
(J) 

X96 

L3 

227 

0.5334 

4671 

2.223 

4.041 

6.999 

X96 

L4 

228 

53.34 

5516 

2.527 

5.004 

3.921 

X96 

L5 

229 

5. 334xl0-J 

4448 

4.318 

1.230 

9.869 

X96 

L6 

230 

0.5334 

5694 

2.667 

4.028 

10.344 

X96 

L7 

231 

53.34 

8452 

5.398 

3.713 

30.834 

34A 

T8 

232 

5.334 

1935 

1.429 

1.541 

1.212 

34A 

T9 

233 

5.334x10 

2157 

5.080 

0.711 

6.799 
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0.943 

9.827 

34A 
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1.611 

11.462 

34A 

T12 
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5. 334x10“** 
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4.763 

0.861 

7.819 

34a 

11 
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1.907 

12.835 

34A 

L2 
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1.706 
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L3 
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1.576 

2.642 
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Til 

247 

5.334x10-2 

2713 
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1.051 

3.632 

584A 

T12 

248 

5.334 

3959 

2.699 

1.681  * 

4.638 

584A 

LI 

249 

5.334x10-4 

2758 

2.159 

1.560 

2.949 
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L2 

250 

5.334x10-2 

3781 

3.493 

1.423 

6.728 

584A 

L3 

251 
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5115 

3.493 

1.842 

9.756 

584A 

L4 

252 

5.334x10-4 

3136 

3.810 

1.101 

6.672 

584A 

L5 
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5.334xl0“2 

4181 

4.445 

1.479 

11.988 

584A 

L6 

254 

5.334 

5560 

5.080 

1.691 

16.965 

584A 

L7 

255 

5.334x10-4 

3247 

8.128 

0.881 

10.214 

88  2A 

T8 

256 

5.334 

■ — 

— 

— 

— 

882A 

T9 

257 

5.334x10-2 

— 

— 

— 

— 

882A 

T10 

258 

5.334x10-4 

2891 

4.445 

1.061 

7.453 

882A 

Til 

259 

5.334x10-2 

3914 

3.175 

1.910 

7.604 

882A 

T12 

260 

5.334x10-4 

3514 

5.080 

0.799 

7.666 

882A 

LI 

261 

5.334 

5783 

5.080 

1.416 

17.982 

882A 

L2 

262 

5.334x10-4 

3781 

3.366 

1.653 

7.700 

882A 

L3 

263 

5.334 

6850 

3.937 

3.121 

16.383 

882A 

L4 

264 

5.334x10-2 

5605 

5.715 

1.615 

19.315 

882A 

L5 

265 

5.334x10-4 

5471 

6.985 

1.217 

22.563 

882A 

L6 

266 

5.334xl0-2 

6939 

5.906 

2.215 

26.269 

882A 

L7 

267 

5.334 

8007 

5.080 

2.802 

25.681 

43 


v 


arxr 


APPENDIX  I (Concluded) 


ENGINEERING  DATA  OBTAINED  FROM  RHESUS 
MONKEY  SINGLE  VERTEBRAL  BODIES 


Specimen 

No. 

VB 

No. 

Displace. 

Rate 

(m/min) 

Ult. 

Load 

(N) 

Def . to 
Ult.  Load 
(mxlO-3) 

Stiffness 

(N/mxl0“) 

Energy  to 
Ult.  Load 
(J) 

W92 

T8 

268 

5.334xl0-4 

867 

2.159 

0.458 

0.859 

W92 

T9 

269 

5. 334xl0-2 

1223 

2.191 

0.631 

1.637 

W92 

T10 

270 

5.334 

1646 

3.651 

0.490 

2.781 

W92 

Til 

271 

5.334xl0-4 

1090 

2.540 

0.494 

1.408 

W92 

T12 

272 

5.334x10-2 

1334 

3.239 

0.480 

2.110 

W92 

LI 

273 

5.334 

1913 

2.858 

1.090 

3.197 

W92 

L2 

274 

5.334xl0-4 

1601 

4.572 

0.457 

4.000 

W92 

L3 

275 

5.334xl0-2 

2535 

4.445 

0.978 

7.355 

W92 

L4 

276 

5.334 

2269 

1.969 

1.168 

1.667 

W92 

L5 

277 

5.334xl0-4 

1668 

1.778 

1.066 

1.347 

W92 

L6 

278 

5.334xl0-2 

1824 

1.905 

1.194 

1.638 

W92 

L7 

279 

5.334 

2402 

1.588 

1.5181.518 

1.678 
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ENGINEERING  DATA  OBTAINED  FROM  RHESUS 
MONKEY  SINGLE  VERTEBRAL  BODIES 
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ENGINEERING  DATA  OBTAINED  FROM  RHESES 
MONKEY  SINGLE  VERTEBRAL  BODIES 
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